Flux period, spin gap, and pairing in the one-dimensional t — J — J'-model 
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Using the factorization of the wavefunction in the t- J- J'-model at smaU exchange couphngs, we 
demonstrate the connection between the existence of a spin gap and an hc/2e flux periodicity of 
the ground state energy. We conjecture that aU spin-gapped SU(2)-invariant Luttinger liquids have 
hc/2e flux periodicity, and that this is connected to the fact that a gapped spin-^ chain always 
breaks translational symmetry by doubling the unit cell. 
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Soon after the discovery of high- Tc superconductivity, 
Anderson proposed that the basic physics of the cuprates 
is that of a doped two-dimensional Mott-insulator H. 
In particular, the Cooper pairs of the superconducting 
state are viewed as the "liberated spin singlet pairs" of 
the insulating host material. While this picture is very 
attractive, it has been difficult to find an explicit model 
for which the proclaimed behavior can be shown to occur 
unequivocally. 

Searching for models of Mott insulators that show su- 
perconductivity upon doping has been the motivation for 
many studies of one-dimensional systems. 0,0, 0,0, IE 
Thanks to methods such as perturbative renormalization 
group and bosonization, considerable knowledge has been 
acquired on the weak coupling phase diagram of both 
strictly one-dimensional _8J and ladder systems .^J. The 
drawback of the weak coupling approach is that it often 
is only an instability analysis. The ultimate statement 
of the quantum phase still rests on certain assumptions 
about the "strong coupling fixed point" of the renormal- 
ization group flow. 

Most strong coupling models cannot be solved analyt- 
ically. A notable exception is the Luther-Emery action 
pH l which describes an electronic liquid with a spin gap 
and dominant singlet-superconducting (SS) correlations 
at large distances. Another interesting analytic method 
for analyzing strong coupling ID models was introduced 
by Ogata and Shiba , and extended in Ref . 3] . This 
method is designed to treat the large U Hubbard model 
(or the small J t-J model). It is based on two facts: 

i) in the limit oiU ^ oo (or J/t ^ 0) the ground state 
of the Hubbard {t-J model) is infinitely degenerate, and 

ii) each of the degenerate states is described by a wave- 
function composed of a product of pure charge and spin 
components For large but finite U (small J/t) one 
can apply degenerate perturbation theory to lift the de- 
generacy. After doing so, the ground state wavefunction 
remains factorized. Moreover, the spin wavefunction is 
given by that of the Heisenberg model on a "squeezed 
lattice" (i.e. the lattice where the unoccupied sites are 
omitted). 

In superconductivity the hallmark of electron pairing is 



the $0/2 = hcjle flux period. In three dimensions, if one 
plots the ground state energy i?($) of a solid supercon- 
ducting torus as a function of the Aharonov-Bohm (AB) 
flux $ through the hole, one finds a periodic function 
with period $o/2. Moreover the energy barrier separat- 
ing the successive minima is extensive. In two and one 
dimensions the flux period is the same. However the en- 
ergy barrier becomes intensive for two dimensions, and 
vanishes as the inverse circumference for one dimension. 

In one dimension the spin and charge degrees of free- 
dom decouple in the low energy and long wavelength 
limit. According to common wisdom, the presence of 
a spin gap implies pairing. It is thus natural to draw 
a connection between the existence of a spin gap and 
a $0/2 flux period. However, since the vector potential 
only enters in the charge action, it is not obvious how 
the presence of a spin gap may affect the flux period. 
The purpose of this paper is to clarify this issue in the 
context of a strongly correlated ID system. 

In the following we study the one-dimensional t-J- J'- 
model, making use the degenerate perturbation approach 
introduced in Ref. 3J ■ The model is defined by the Hamil- 
tonian 



clcrCi+l.cr + h.C. )V 
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(1) 



describing N electrons on a ring of L sites in the presence 
of an Aharonov-Bohm flux $. Here, the projection oper- 
ator V excludes states with doubly occupied sites, and 
the Si are spin-1/2 operators. At J = J' = the model 
corresponds to the U — 00 Hubbard model. As pointed 
out in Ref. > in this limit the eigenstates factorize into 
products of pure charge and spin states. This property 
has been used extensively to study the large (but finite) U 
Hubbard model [ill |l3l| , which is related to the small 
J t-J model. Much less analytic work has been done on 
the t-J-J' model with a finite a = J' /J, because the 
model is no longer integrable. In this case, however, the 
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degenerate perturbation approach introduced in Ref.|3| 
still allows one to determine the ground state properties. 

In the following we will use this method to study the 
ground state energy of J^l as a function of the AB-flux 
$. We begin by defining the TV-particle wavefunction of 
the system 

l*> (2) 

on the domain 



D 



. . . xjv) e xi<X2< - ■ ■Xn<xi+ L} . 



Here j^*) and |0) are the state of the system and the 
vacuum of the fermionic operators Cx^a- = Cx+l.o-, re- 
spectively. The fermion antisymmetry and the periodic 
boundary condition imply 

5' (xi,(Ti, ...,XN,aN) = 

{-l)(^~^)^Sj{x2,a2,...,XN,<JN,ixi+L),ai) (3) 

At J = J' = 0, each eigen wavefunction of |^ factor- 
izes into a product of a charge and a spin wavefunction 

ini, 

^{Xi,ai, ...,XNj(Tn) = f{xi . ..XN)g{cri ■ . .CrAr)(4) 



where 



f{xi . ..xn) 



det [exp(ifciXj) ] 



(5) 



is a Slater determinant constructed from TV plane waves. 
It is of central importance here to observe that for a 
finite ring with periodic boundary conditions, the spin 
part and the charge part of the wavefunction are not 
completely independent. Specifically if we quantize the 
N momenta in Q according to 



2tt K 
where qj G Z and K £ [0, 27r), 



(6) 



then the condition (j2Jl requires the spin wavefunction to 
satisfy 

g(cri . . . o-at) = e'-^ g((T2 . . . CTAT, cri), (7) 

which implies that K is the "spin momentum" on the 
squeezed lattice. Thus by means of the momen- 
tum K of the spin wavefunction injects a twist into the 
charge wave function. For large but finite U the exact 
ground state wavefunction of the Hubbard model remains 
of the form given by Eq. Q), and the same relation be- 
tween charge twist and spin momentum is observed 
Within the degenerate perturbation approach introduced 
in Ref.Q the same still applies to the ground state of Q 
for any value of a = J' / J in the limit of vanishing ex- 
change couplings. In this limit, all solutions of the form 
Eq. are degenerate in the spin wavefunction, which is 
only required to have the spin momentum K determined 
by the twist of the charge wavefunction. 



To first order in the exchange couplings, this degen- 
eracy is lifted by an effective Hamiltonian acting in the 
"squeezed" space of N spins 0: 



Heff = ]^ Jeff ^ Sj ■ Sj + i + Jiff ^ Sj ■ Sj+2 

j=i i=i 



Jeff = J / + J' - nj+i)ni+2) f 

J'^ff = J' {nini+ini+2) f (8) 

Here, {) f denotes a spinless fermion expectation value 
with respect to the wavefunction / displayed in 10}. 

We now focus on the case of constant a where a > 
dc ~ -241 (14] . In this regime numerical and analytical 
works suggest the phase diagram shown in Fig. ^ At 
zero doping the spin chain corresponding to the model 
at half filling (TV = L) is gapped. The spin gap will 
survive for a finite range of doping a; = 1 — N/L < Xc 
(Fig. nj, and the effective spin Hamiltonian © may be 
used to calculate Xc exactly in the limit J/t —t 0. This 
was first proposed by Ogata et al. 0] and has been con- 
firmed numerically in Ref. [isj . Thus the validity of 
degenerate perturbation theory is well established. (We 
will give a detailed discussion of the involved subtleties 
elsewhere (see also j?!)). 

To first order in J/t and J'/i, the ground state energy 
of Eq. takes the form 



Etnt — Er 



E, 



(9) 



where Ec is the kinetic energy associated with the charge 
wavefunction /o, and Eg is the ground state energy of © 
with the spin momentum K. For any given AB-flux, /o 
will be of the form given in 101, where the TV consecutively 
occupied momenta kj are given by © with 

1, J-1...TV (10) 



qj = To + J 
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FIG. 1: Sketch of the zero temperature phase diagram of 
as obtained in 0, ITsl for q = i . The spin-gapped region is 
divided by a crossover (dashed) between regions of dominant 
singlet superconducting (SS) and charge-density-wave (CDW) 
correlations. 
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FIG. 2: Ground state energy according to Q, Ulljl for x < Xc and x > Xc- The energy scale is t/L. a.) x < Xc (spin-gapped 
case). Periodic pattern of branches with K = and K = tv separated by half a flux quantum. The thick line represents the 
ground state energy at given $. The ground state momentum P changes by kf between adjacent branches, b) j: > a;c (no spin 
gap). The i$o-periodicity is destroyed by a relative shift of order J/L between K = and K = tv branches. 



and qo is an integer. The kinetic energy is given by 
E,{qo,K,<f) = -2^^ cos (kj + 

sin(fcf) /, 27r K - n 27r $ \ 

= -2t V-^ cos kf H On H ! 

sin(f) V L L L $0/ 

(11) 

In Eq. ifTT)) kf is defined as fc/ = ttN/L. Note that for all 
charge wavefunctions characterized by Eq. ((BJ, (|10|l with 
different go and K, the effective couplings appearing in 
© are the same. 

We will consider an even number of particles N from 
now on. In the spin-gapped regime < x < Xc, the 
effective Hamiltonian (0 has two lowest-energy states 
with spin momenta K — and K — tt. These two states 
are separated from other spin states by an energy gap of 
order J or J'. For a finite ring, the energy difference be- 
tween these two lowest spin states vanishes exponentially 
with the circumference of the ring. Thus in the thermo- 
dynamic limit these two states become degenerate. 

For fixed we choose K and qo so that the total 
energy is minimized. This minimization can be achieved 
by first minimizing Ec by varying qo for fixed then 
minimizing + Ec with respect to K. When the first 
minimization is achieved the argument of the cosine is 
always of order 1/L regardless of the values of K. Hence 
the effect of varying K in the second minimization can 
only result in 0{t/L) modulations in Ec- We consider 
the limit t/L <^ J here. It then follows that the value 
of K must be either or tt. Other choices of K would 
increase Es by the spin gap of order J, which cannot be 
compensated by the possible lowering of Ec- 

After substituting the optimum value of go for i^T = 
or = TT back into Eq. we obtain two branches of 
energy versus $ curves shown in Fig. \^). The lower 
envelope of these curves is the ground state energy as 
a function of $ for small J/t and J' /t. One observes 
that this function does indeed show a period of $o/2, 
owing to the existence of two different types of branches 
corresponding to if = and K = t:, respectively. This 
structure is resemblant of that proposed for the dimer 



model jlTl Il8l |. It is also interestirig_to note that K — n 
is the analogy of the "vison" flux 19] in ID. 

The situation is fundamentally different in the regime 
X > Xc (Fig. P), where the spin gap vanishes. Here the 
ground state of the effective spin Hamiltonian ^ has 
K = Q or K = TT depending on whether N = 4m or 
N = Am + 2 (see e.g. "20]). For a finite chain of length 
N = 0{L) the first spin excited states at A' = vr or 
K = {i. e. whose momenta differ by tt from that of the 
respective ground state) have excitation energies of order 
J/L. This energy gives rise to the relative shift between 
the K = and K — tt branches shown in Fig. Eb) 
for N — 4m. The resulting lower envelope is illustrated 
for a small but finite J/t, where the flux period is now 
$0- For iV = 4m -I- 2 the shift between the K = and 
K = TT branches is opposite in sign. This behavior is well 
demonstrated in the large repulsive U Hubbard model 
|2l| . where no spin gap is present. 

It is sometimes felt that the existence of metastable 
minima of £'($) at $o/2 intervals is the sign of a pairing 
tendency, even though $o/2 is strictly not the flux period. 
FigEb) presents a clear counter example of this type of 
reasoning. Indeed according to Fig. \3p) this can hap- 
pen in a regime of the phase diagram where the ground 
state is neither spin-gapped, nor features a dominance of 
superconducting pairing correlations. 

On the other hand, for x < Xc, there is a spin gap but 
no charge gap, and the ground state energy is a periodic 
function of $ with period <I'o/2. It has thus all the char- 
acteristics of a superconductor. However for J/t ^ 1 
the superconducting correlations (SS) are weaker than 
the charge density wave (CDW) correlations (Fig. 
In this regime we can think of the system as being close 
to a superconductor-(Cooper pair) insulator transition 
due to strong quantum fluctuation of the phase of the 
superconducting order parameter. Here, a weak exter- 
nal perturbation such as disorder can easily localize the 
Cooper pairs and drive the system insulating. As the 
system crosses the crossover line in (Fig. the phase 
fluctuations become much less severe so that the SS cor- 
relations become dominant over the CDW correlations. 
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We expect the appearance of a $0/2 flux period to hold 
in the entire spin gapped regime. Indeed, more detailed 
considerations show that the arguments given here are 
not limited to first order perturbation theory |l6|. In 
particular, we find that second order contributions to 
the ground state energy may be incorporated into the 
efl'ective spin Hamiltonian. The latter will then also de- 
pend on the twist of the zeroth order wave function and 
on flux, yet Heff{K, $) = Heff{K - tt, $ + $o/2) con- 
tinues to hold. Furthermore, we have shown that all 
the results presented here also follow from a weak cou- 
pling/bosonization procedure |2^. 

The analysis presented here can be applied to a wide 
class of models of the form Eq. where the second line 
is replaced by a more general spin-chain type of Hamilto- 
nian. Most features of the phase diagram shown in Fig. 
^will likely survive as long as the spin chain at half fill- 
ing is gapped. In particular, the spin gap will survive 
for a range of doping, and phase separation will occur 
at sufficiently large values of J/t, where J is an appro- 
priate energy scale for the spin couplings. As the phase 
separation line is approached, the charge compressibil- 
ity diverges, and Luttinger liquid physics then implies a 
regime of dominant SS correlations. Our analysis on 
flux period will then carry over to this more generic case, 
provided that the gapped spin state at a; = also breaks 
translational symmetry by doubling of the unit cell, anal- 
ogous to the dimerization that occurs in the J- J' model 
at half flUing. Such an example is given by the t-Jz model 
studied in Ref. [13 ■ Although in Ref. the possibil- 
ity of ^o/n flux periods {n > 2) has been postulated for 
models of the type considered here, only the case n — 2 
has been found for the t-J^ model. This follows easily 
along the line of arguments given here, and we believe 
that only n = 1 and n = 2 are found in generic models. 

If, on the other hand, a gapped spin-i chain exists that 
does not break translational symmetry, it appears that a 
doped model with a spin gap could be constructed which 
does not feature $o/2 flux quantization as displayed in 
Fig. However, such a state would violate the Lieb- 
Schultz-Mattis theorem (23|. For SU(2)-invariant spin- 
i chains in one dimension, we are aware of only one 
way to create a spin gap, i.e. breaking the translational 
symmetry by doubling the unit cell. Hence there seems to 
be an intimate relation between this fact and the possible 
universality of the <E>o/2 flux period which we postulate 
below. 

We note that a $0/2 flux period associated with a spin 
gap has also been observed in numerical studies of a two- 
leg ladder [i^. This suggests that our main conclusion 
may be generalized beyond the purely one-dimensional 
case. However, a two-leg ladder has an even number of 
sites per unit cell. Here the undoped system may have 
a spin gap due to the formation of singlet pairs located 
on the rungs, which does not require symmetry breaking. 
These singlet pairs become mobile upon doping, and the 



above notion of symmetry breaking in some internal spin 
space is not required to explain the $o/2 period. 

To conclude, we have demonstrated the relation be- 
tween a $o/2-periodicity in the ground state energy and 
the existence of a spin gap in the small exchange limit 
of the t-J-J' model. Based on these flndings, we con- 
jecture that the observed $0/2 flux period is a universal 
property of spin-gapped SU(2) invariant one-dimensional 
systems of spin- ^ particles with gapless charge degrees of 
freedom. In particular, the value of the charge Luttinger 
parameter is not a determining factor of the flux period- 
icity, as our result did not require the predominance of 
singlet superconducting correlations. Our findings fur- 
ther suggest an intimate relation between the proposed 
universality of the $o/2 flux period and the fact that 
all gapped SU(2)-invariant spin-i chains feature broken 
translational symmetry with a doubling of the unit cell. 

This work has been supported by DOE grant DE- 
AC03-76SF00098. 
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